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1  INTRODUCTION 

This  study  is  intended  to  facilitate  a  gravity  field  adjustment  containing 
observations  of  geo  id  undulations  along  profiles.  The  resolution  of  the 
adjustment  is  assumed  to  he  represented  by  an  (n.n)  spherical  harmonic  expansion 
or  its  equivalent  This  implies  that  the  geoidal  detail  heyond  the  (n.n)  model 


enters  the  realm  of  observational  "errors' 


It  acts  as  if  increasing  the 


variance  of  the  observational  noise  proper  and  Introducing  (or  increasing) 
invariances  Such  an  added  effect  can  be  evaluated  with  the  aid  of  the  geoida) 
covariance  function,  and  it  leads  to  the  final  variance  covariance  matrix 
associated  with  geoidal  observations  In  a  straightforward  approach  this  matrix 
would  be  inverted,  resulting  In  a  fully  populated  weight  matrix. 

At  first  sight  ,  one  might  object  to  "'mixing'’  quantities  such  as  the 
observational  noise  proper  and  the  neglected  geoidal  detail,  where  the  latter 
lias  the  nature  of  a  (neglected)  signal.  However,  their  simultaneous  treatment 
is  compatible  with  the  collocation  principles  pioneered  by  Moritz  [1980]  and 
used  here  in  forming  the  above  geoidal  variance  covariance  matrix.  Under  the 
collocation  operator  MsEM;ME,  where  the  operator  M  symbolizes  the  average  over 
the  unit  sphere  (applied  to  the  signal)  and  the  operator  E  symbolizes  the 
mathematical  expectation  (applied  to  the  observational  noise  and  other  random 
effects),  the  distinction  between  the  original  characteristics  of  such 
quantities  is  essentially  obliterated. 

A  typical  example  of  geoidal  observations  along  profiles  is  offered  by 
satellite  altimetry.  In  considering  the  thousands  of  orbital  arcs  (profiles) 
entering  a  global  adjustment,  a  full-scale  inversion  of  the  variance  covariance 
matrices  associated  with  these  arcs  would  be  computationally  prohibitive.  On 
the  other  hand,  the  neglect  of  all  covariances  would  be  unacceptable  from  the 
theoretical  standpoint  A  compromise  solution  was  offered  in  Chapter  4  of 
(81aha,  1984],  which  resulted  In  a  tri  diagonal  weight  matrix.  In  that 
approach,  the  weight  matrix  was  formed  almost  as  efficiently  as  in  t  e  diagonal 
case  (with  neglected  covariances),  while  the  corresponding  errors  in  the 
covariance  function  were  lowered  to  only  about  40%  of  those  characterizing  the 
diagonal  case  The  standard  of  comparison  was  the  case  called  rigorous,  where 
the  covariances  are  properly  accounted  for  and  the  weight  matrix  is  fully 
populated  The  present  analysis  is  guided  by  the  need  to  reduce  such  errors 
even  more  significantly 


1 


The  basic-  requirement  leading  to  highly  patterned  characteristics  of  the 
variance  covariance  and  the  weight  matrices  is  a  uniform  spacing  of  observations 
along  the  profiles  With  regard  to  altlmetrlo  data,  this  requirement  presents 
no  difficulty  Although  the  effect  of  the  observational  noise  proper  and  other 
sources  of  error  (eg.  unmodeled  tidal  variations)  was  estimated  in  [Biaha, 
1984)  and  added  to  the  variance  computed  from  the  geoidal  covariance  function, 
here  this  effect  will  be  disregarded.  For,  as  is  evident  from  this  reference, 
not  only  is  its  estimation  very  approximate,  but  it  modifies  the  geoidal 
variance  by  only  a  relatively  insignificant  amount.  However,  if  desired,  this 
effect  and  other  modifications  such  as  scaling  of  the  covariance  function 
can  be  included  without  changing  anything  in  the  analysis.  Likewise,  the  same 
analysis  is  valid  for  any  (n,n)  model  describing  the  adjustable  gravity  field, 
for  any  length  of  individual  geoidal  profiles,  and  for  any  uniform  observational 
i nterva 1 . 


«*  w-*  v<v 


2.  GEO I DAL  COVARIANCE  FUNCTION  AND  ITS  MODIFICATIONS 

The  formation  of  the  weight  matrix  of  observations  will  be  carried  out  in 
lour  cases  The  first  case  is  represented  by  the  unaltered  covariance  function. 
The  second  case  is  represented  by  a  modified  covariance  function  as  obtained 
through  a  fit  by  two  parameters,  and  a .  The  third  case  is  characterized  by 
a  more  profoundly  modified  covariance  function  as  obtained  through  a  fit  by  one 
parameter,  a^.  And  in  the  fourth  case,  all  covariances  are  neglected  while  the 
variance,  computed  from  the  geoidal  covariance  function,  is  kept  intact. 

The  representation  of  the  covariance  function  for  geoid  undulations  is 
based  on  t fie  results  of  [Tscherning  and  Rapp,  1974] .  These  results  were  adopted 
in  [Klaha,  1982]  to  give  the  k-th  degree  variance  (in  meters  square)  as  follows: 

=  17981x0. 999617k+2/[(k-l)(k-2)(k+24)]  . 

The  covariance  function  for  geoid  undulations  corresponding  to  the  truncation  of 
the  spherical -harmonic  model  at  the  degree  and  order  (n.n)  is  then  given  as 

0 (<l>)  =  M(N,N'  }  =  l  a *  Pk(cos *)  , 

where  the  summation  is  carried  out  from  the  degree  n*l  to  infinity.  The 
operator  M  indicates  averaging  over  the  unit  sphere,  N  and  N'  are  the  geoid 
undulations  at  any  two  points  P  and  P'  separated  by  the  spherical  distance  ifr , 
and  Pk(cosii)  are  the  Legendre  polynomials  in  the  argument  cos 0.  Since  Pk (1)^1. 
one  has  for  the  geoidal  variance,  D(0),  imputable  to  the  neglected  degrees  (from 
n + 1  to  infinity) : 

0(0)  =  l  ak  . 

In  practical  computations,  the  degree  "infinity"  is  substituted  for  by  1000  or 
some  other  sufficiently  large  number. 

The  truncated  model  in  this  analysis  is  represented  by  (22,22)  spherical 
harmonic  expansion,  i.e.,  n=22.  With  this  truncation,  the  geoidal  variance  is 

a 2  i  0(0)  --  11 .367  m2  .  (1  ) 

The  observational  interval  is  stipulated  to  be  0.5°.  The  geoidal  covariances, 

2 

including  the  variance,  are  presented  below  in  a  direct  form  (the  units  are  m  ), 

2 

as  well  as  in  multiples  of  a  ,  in  which  case  the  entries  are  denoted  "p" . 
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<l>  : 

0 

O 

tr> 

o 

1 .0° 

1.5“ 

2.0° 

2.5° 

3.0° 

l)  ; 

1  1 . 367 

10.661 

9.293 

7.693 

6.047 

4.462 

3.003 

P  '■ 

1 

.93789 

.81754 

.67678 

.53198 

.39254 

.26419 

For  the  distances  ip^3.5°  through  iji=80  the  values  of  p  would  be  .15035,  .05296, 
.02710,  .09000,  .13627,  .16715,  -.18413,  -.18888,  -.18316,  -.16917. 


The  second  case  is  based  on  the  Model 

Pk  -  “iVl  +  a2Pk-2  ' 

where  k^l  and  where,  by  definition. 

p0  =  1  '  p(  1)  =  P1  ' 

The  value  pQ=  1  in  (3b)  corresponds  to  p=  1  for  t^-0  in  (2).  Equations  (3a-c) 
agree  witli  the  Yule  scheme  as  is  apparent,  for  example,  from  (6.21)  in  (Kendall, 
1973],  The  third  case  is  characterized  by  a2=0,  f°urth  case  is 

characterized  by  a^a^O. 

From  (3a~c),  the  values  of  p  for  k>l  can  be  expressed  in  terms  of  the 
parameters  a  and  a as 

2 

p\  =  “l^1  "a2*  ’  p2  =  a2  +  ^/(1-a^)  - 
P3  -  (a,/(l-«2)JC«f  +  2a2  '  a2]  . 

which  is  valid  for  all  three  cases  (second,  third,  and  fourth).  A  better 
formulation  for  the  second  case  is  given  by 

Pk  =  pksin  (k0  +  (3)  / sin$  , 

valid  for  k?0,  where 

p  -  (  “a2^1/2  ’  cos0  =  a1/'2P  •  tg/3  =  [(1 -o2)  tgfll/(l +  a2)  . 

This  formulation  can  be  derived  using  an  approach  similar  to  that  on  page  73  of 
[Kendall,  1973],  From  trigonometric  considerations  it  then  follows  that 

0  <  -«2  $  1  ,  0  (  jj  <  2(-<*2)1/2  <  2  .  (4a. b) 

Equations  (4a, b)  serve  to  verify  the  admissibi  1  lty  of  and  c*2 


(3a) 


(3b,  c) 


I 


Upon  using  (3a)  and  a  similar  equation  for  p 


t he  solution  of  these  Iwi 


k-  1 

equations  in  two  unknowns  (a^  and  a^)  yields  for  the  second  cast 


'Vi'k 


V2Vll/d  • 


°2  ‘  ''k  l'k.l 


4*'*  ■ 


d  =  p 


k  1 


Pk-2Pk  ’ 


where  the  values  for  pQ  and  p  ^ 
has 


( r.a  ,  t. ) 
(5<  ) 

are  given  by  (3b, c).  For  the  third  case,  one 


«i  *  Wi  •  “2 ‘ 0  ■ 

where  (3b)  remains  valid,  while  the  fourth  case  has  already  been  described  by 
-  0  ,  =  0  .  (7a, b) 

In  evaluating  the  second-case  results  and  a 2  for  a  given  profile,  one 
could  use  all  of  the  available  p's  in  <5a-c)  and  average  the  results.  However, 
here  such  averages  are  used  only  in  the  role  of  initial  parametric  values  in  a 
least-squares  adjustment.  Upon  considering  n*l=17  values  of  p  associated  with 

ip-0  through  ii>=8°  in  (2)  and  below  it,  denoted  pQ.  pJ(  p^ .  p^,  it  follows 

from  (5a-c): 

a1  =  1.70163  ,  «2  =  -.75674  ,  (8a, b) 

where  n=16  values  of  and  have  been  used  to  compute  either  average.  (The 
symbol  n  need  not  be  confused  with  the  same  letter  which  earlier  indicated  the 
degree  and  order  of  the  spherical-harmonic  model.)  With  the  values  in  (8a, b), 
equations  (4a, b)  are  satisfied.  In  the  third  case,  the  initial  value  of  has 
been  computed  by  averaging  only  the  first  six  values  in  (6a, b),  resulting  in 

«1  =  .80573  a2  =  0  .  (9a,b) 

The  reason  for  this  limitation  stems  from  the  fact  that  beyond  tl>= 3°  the  values 
of  p  become  small  and  even  change  the  sign. 

The  adjustment  model  in  both  the  second  and  the  third  cases  is  general,  in 
the  sense  that 

F(Xa, L3)  =  0  , 

3  3 

where  X  symbolizes  the  adjusted  parameters,  and  eventually  a 2>  and  I, 
symbolizes  the  n  adjusted  "observations”,  p^ ,  p2 .  Pn .  Since  pQ-l  is  a 
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ti 


with  I  l.‘  t  tie  ntisei  ved  values  of  p  represented  tieie  ti\  t  hr  rigorous  ,  ,i  n  <■ 
t  J  )  Itn!  lie  i  111*  it.  With  t  tie  exception  of  the  first  Ip  !  ,  !  fit-  i  et  t  v  ,  ten  '  >■  ■ 

!  t 

f  he  id  jus  t  men  t  rot  m  t  ions  to  t  tie  Initial  parametric  mines  >  and  the  .  r  • 

X'  denotes  the  residuals  The  least  squares  solution  i  ri  this  getter  a  I  <*,1  ms  r  men* 
me t  tu it  is 


(  A  F  M  '  A  )  'aTM  '  W  . 


i  io 


where 


M  BIB 


!  i 


L  being  t  tie  variance  covariance  matrix  of  the  observations  .  taken  here  as  I 


the  second  case,  the  general  node!  follows  from  ('la)  as 


V'k 


<Vk  2  =  ° 


l  12  ! 


wliei  e  t  tie  superscript  "a"  has  been  omit*ed  For  k*3,  the  k  t  h  row  l  ri  t  tie  matrix 
A,  the  three  nonzero  elements  in  the  k  th  row  of  the  matrix  B,  and  the  k  t  fi 
element  In  the  vector  W  are,  respectively, 


°2  “l  1] 


pap  a  o 
Mk  1 pk  1  2'k  2 


(  1  t<i  ,  ft  .  i  \ 


‘  r  k  1  f  k  2‘  ' 

where  the  superscript  "b"  in  conjunction  with  the  p's  has  also  been  omitted 
The  element  in  (13b)  is  located  at  the  diagonal  of  B  For  k  1.  equation 

(13a  r)  become 


[  1 
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[1  «2  0  0]  , 


V1  V  a!  ' 


<13  a  h  i 


whereas  for  k  2  they  become 


[  b,  1 ]  . 


(  «j  1  0]  , 


o  a  o  <x 
ui°i  u2 


(  1  3 '  a  ,  b ,  r  ) 


•'  6  and  I*  ar e  seen  in  be  (  n*  2  )  ,  (  n  *  n  )  ,  and  In-  1  ) 

1  b’  ’  b  i  l  <1  i  isi-  (.in  be  described  in  A  pel  fee  t  analogy  to  tin 

bi’  ’ se:  oini  i  u  I  uin  in  A  is  nonexistent  .md  that  a  0 

2 

®l"  .  n«:  *  be  so  (nr  inn  X  in  (10).  one  can  t  nke  advantage  of  tin-  (net 
t'!!':  '.  M  '  •in  be  inverted  analytically  Since  H  is  a  regular 
•  c -  00!  nnns l ngu lari,  it  follows  from  (11)  that 


sim;  :  i  f  red  by  l  I  Due  to  the  patterned  characteristics  of  B 

'  1  fs't  M  t  b).  it  is  readily  verified  that  the  matrix  B  1 

i  **•  :  •  :  ’  ■! .  r ,  >;  t..  'he  following  algorithm  Initially,  a  "dummy"  column 

if*  '  w r  •  '  -  i  element  is  d  1.  second  element  is  d  a  and  any 

!  ■-  1 

. me  u  *  .  s  d  -id  -ad  1-3,4 . n  Next,  the  first  column  in 

!  1  1  2  1  3 

'  n-  '  ■  t,e  dummy  loiuam  times  the  factor  1  ( ]  -a  )  .  Finally. 

i  *  i.»  i :  s  g  i  ■  n  a !  the  second  column  in  B  1  consists  of  the  first  n-  1 

•  • ■  ""  1  • he  lummy  ■  o i  iimn  Similarly,  starting  at  the  diagonal,  the  third 

.■a*-  i  u  i  s  .  s  r  s  of  the  first  (n  2)  values  of  the  dummy  column,  and  so  on, 

'  ’b"  -ii agonal  element  1  In  the  last  column  is  filled  in.  The  elements 

>b<.vi  t  !.i  d  •;  ur'iti.i  i  in  these  columns  are  zero,  i .  e .  .  the  matrix  B  1  is  lower - 
ii'K’ilu  The  1 1,  ]  r  d  case  again  follows  by  setting  a2  =  0. 

\  nearly  perfert  i onvergence  in  the  iteratec  adjustment  (10)  has  been 
achieved  in  fom  iterations  in  both  the  second  and  the  third  cases.  The  results 

*  '  '  !c  sei  ond  case  ire 

a,  1  69244  ,  a2  76172  ,  (14a, b) 

satisfying  (4a  hi  The  i  17  values  of  the  adjusted  "observations"  p 
■or  i  espond  i  eg  to  !2)  arid  its  extension  are 

e'sei.uidi  1  96066 .  86417,  73080,  .57858,  ,42255,  .27442, 

14258.  03228.  05397,  11594,  .15510,  .17419,  (15a) 

1  ~l,(,7  .  166  n  ,  14690,  .  12194  . 

1  nmpared  to  (2)  anti  its  extension,  the  root  mean  square  of  the  errors  in 

I)1  ser  ond  )  1  s 

r ms ( second )  02989  (15b) 
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Similarly,  the  third  case  yields 


=  .80243  ,  =  0  , 

and 


(16a, b) 


p  (third)  =  1,  .80243,  .64389,  .51668,  .41460,  .33268,  .26695, 

.21421,  .17189,  .13793,  .11068.  .08881,  .07126,  .05719, 

.04589,  .03682,  .02955  ; 

rms (third)  =  .16907  . 

The  fourth  case  is  trivial,  namely 

p(fourth)  =  1,  0,  0,  0 . 0  ; 

rms (fourth)  =  .40040  . 

In  terms  of  the  root-mean-square,  the  error  in  the  third  case  amounts  to 
approximately  42*  of  the  fourth-case  error,  in  line  with  the  result  achieved  in 
fBIaha,  1984]  for  the  (14,14)  truncation.  However,  the  error  in  the  second  case 
is  much  smaller,  amounting  to  under  8*  of  the  fourth-case  error.  The  second- 
case  error  is  only  about  18*  of  the  third-case  error.  Such  a  significant 
reduction  indicates  that  no  further  refinements  are  likely  to  be  necessary  in 
the  future.  Indeed,  the  covariance  function  itself  is  only  approximate,  and, 
furthermore,  the  values  it  produces  are  indicative  of  a  global  situation,  not 
that  pertinent  to  a  particular  geoidal  profile. 

On  the  other  hand,  it  may  be  beneficial  to  repeat  the  above  adjustment  for 
varying  lengths  of  the  profile,  as  well  as  for  a  different  observational 
interval.  For  example,  the  current  results  (14a, b)  may  be  deemed  satisfactory 
for  a  profile  length  of  6°-10°,  whereas  for  a  different  length  the  final  values 
of  and  may  slightly  change.  But  the  parameters  adjusted  by  the  general 
method  for  any  length  of  the  profile  are  preferred  to  the  approximate  values 
obtained  by  averaging.  In  the  situation  under  scrutiny,  the  approximate  values 
in  (8a, b)  would  produce  an  rms  error  over  three  times  as  large  as  that  presented 
in  (15b),  obtained  with  the  adjusted  values  (14a, b). 


(17a) 


(17b) 


(18a) 
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3.  HEIGHT  MATRICES 


If  the  coefficients  p^  fulfilled  (3a-c)  exactly,  the  matrix  (19)  could  be 
inverted  analytically,  resulting  in  the  weight  matrix  P  for  the  second  case 
presented  below,  where  the  factor  q  is  given  by 


va 

n 


P  =  (1/ff  )q  I  0 


1 

~ai 

®2 

0 

0 

0 

'al 

l.a2 

1 

-otjd-ag) 

"a2 

0 

0 

*2 

-<*1  ( l-a2) 

,  2  2 
1+Va2 

^l(l-a2) 

-<X2 

0 

0 

“a2 

-ai(l-a2) 

l+a2+a2 

1  2 

l 

» 

1 

ft 

CO 

‘a2 

(21b) 

The  entries  In  all  the  rows  except  the  first  two  and  the  last  two  repeat  in  the 
fashion  Indicated.  The  last  two  rows  are  symmetric  with  respect  to  the  first 
two.  The  Matrix  P  in  the  third  case  would  be  formed  exactly  as  in  (21a,b), 

except  that  a  would  be  replaced  by  zero.  And  the  trivial  fourth  case  would 

*  2 
have  both  and  o 2  replaced  by  zero,  resulting  in  P=(l la  )I. 

In  the  practical  procedure  suggested  herein  the  coefficients  p  given  by  the 
geoidal  covariance  function  are  modified  to  conform  to  (3a-c).  As  (21b) 
indicates,  the  weight  matrix  P  is  then  five-diagonal,  all  the  other  entries 
being  zero.  With  «2=0  this  *atrix  would  be  tri-diagonal,  and  with  a  =a  “0  it 
would  be  diagonal.  Accordingly,  the  first  case  will  also  be  called  "rigorous", 
while  the  second,  third,  and  fourth  cases  will  also  be  called  "five-diagonal", 
"tri  diagonal",  and  "diagonal",  respectively.  In  [Brown  and  Trotter,  1969],  the 
tri  diagonal  and  the  f i ve -d iagonal  cases  were  referred  to  as  the  first-order  and 
the  second  order  autoregressive  processes,  respectively.  It  is  noted  that  for 
the  purpose  of  that  reference,  the  factor  q  was  omitted.  Although  the  present 
approach  could  be  extended  to  a  seven-diagonal  case,  etc.,  the  f ive -diagonal 
weight  matrix  is  deemed  completely  satisfactory  for  geoidal  profiles.  According 
to  the  discussion  at  the  close  of  Chapter  2,  the  errors  committed  by  modifying 
the  covariance  function  in  this  case  are  negligible. 


*  -  * 


The  validity  of  (21a, b)  can  be  confirmed  upon  performing 


PZ  -  I  ,  (22) 

where  the  entries  pk  in  (19)  are  taken  from  (3a-c).  In  considering  the  first 
row  of  P,  the  first  row  of  (22)  follows  from  the  following  three  identities: 


q(l -«1i°1-«2P2)  =  1  ■  (l-a2)p][  -  ftj  =  0 


Pi  +  2  alPi  +  l  a2pi  °  ’ 


where  i>0.  The  second  and  the  third  identities  correspond  to  (3a-c).  The  other 
elements  in  (22)  could  be  verified  in  a  similar  fashion.  For  example,  the 
element  (2,2)  is  the  consequence  of 

and  further  diagonal  elements  (except  for  the  last  two)  are  the  consequence  of 

q[”a2^2ai^  a2)pl+^1+ai+0I2^_ai^  a2^Pl  0I2P2^  =  1  ' 

The  same  confirmation  is  valid  also  with  a  =0  and  with  a  =a  =0. 

w  1  Ct 


4.  SIMPLE  NUMERICAL  EXAMPLE 


In  this  example,  m  =  17  observations  of  geoid  undulations  are  considered 
distributed  in  0.5°-intervals  spanning  a  8°-profile.  From  (2)  and  its 
extension,  one  forms  the  rigorous  (17x17)  variance-covariance  matrix  as 


r 

1 

.93789 

.81754  .  . 

.  -.16917 

.93789 

1 

.93789  .  . 

.  -.18316 

I  =  11.367 

.81754 

.93789 

1  .  . 

.  -.18888 

i 

-.16917 

L 

-.18316 

-.18888  .  . 

1 

The  inversion  of  this  matrix  yields  a  fully  populated  weight  matrix  P,  computed 
with  the  aid  of  the  Choleski  algorithm  but  not  listed  here. 

The  variance-covariance  matrix  modified  according  to  the  second-case 
principle  would  follow  from  (15a).  However,  it  is  not  needed  explicitly  since 
the  five-diagonal  weight  matrix  Is  given  directly  by  (21a, b)  with  the  numerical 
values  of  and  a 2  presented  in  (14a, b).  Similarly,  the  tri-diagonal  weight 
matrix  of  the  third  case  follows  from  (21a, b)  with  and  (*2=0  given  by  (16a, b). 
Finally,  the  diagonal  weight  matrix  of  the  fourth  case  is  P  =  ( 1/11 .367)1 . 

The  17  observations  considered  are  used  to  adjust  one  parameter,  taken  here 
as  a  point-mass  magnitude.  The  point  mass  itself  is  chosen  to  have  a  central 
location  with  respect  to  the  profile,  and  to  be  located  at  the  depth  d=704km 
beneath  the  earth's  surface  represented  by  a  sphere  of  radius  R=6371km.  Thus, 
the  central  angle  ^  varies  between  0  and  4°  on  either  side  of  the  point  mass. 

The  mathematical  model  for  geoid  undulations  is  written  as 

Ni  =  (1/G) (kM) j/sij  .  1=1,2 . m  , 

where  G  is  the  average  value  of  gravity  adopted  as  980gal;  (kM)j£x  is  the 
unknown  parameter,  i.e.,  the  gravitational  constant  times  the  mass  of  the  point 
mass  "j";  arid  Sjj  is  the  distance  between  the  observational  point  "i"  and  the 
point  mass  "j".  This  distance  is  computed  as 
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where  ip  ^  is  the  angle  ip  corresponding  to  the  above  points  1  and  j. 

The  17  observations  of  geold  undulations  (In  meters)  are  simulated  In  a 
symmetrical  manner,  the  first  and  the  last  being  zero  and  the  middle  one 
directly  above  the  point  mass  being  5m,  as  follows: 

0  000,  0.975,  1.913,  2.778,  3.536,  4.157,  4.619,  4.904,  5.000, 

4.904,  4.619,  4.157,  3.536,  2.778,  1.913,  0.975,  0.000  . 

These  values  represent  the  elements  of  the  vector  called  e.  The  units  are 

_  q 

chosen  such  that  x  is  expressed  in  terms  of  10  times  the  earth's  kM  (the 
latter  is  3 . 986005x io14m^/s^ ) .  The  design  matrix  A  has  the  dimensions  (17xi). 
Its  elements,  symmetric  about  the  middle,  are 

.49635,  .51232,  .52748,  .54143,  .55370,  .56384,  .57144,  .57615,  .57775, 

.57615,  .57144,  .56384,  .55370,  .54143,  .52748,  .51232,  .49635  . 

The  parametric  least-squares  adjustment  yields 

x  =  (ATPA)  1ATPe  , 

a2  =  (ATPA)  1  ; 
x 

the  symbol  x  denotes  here  the  least-squares  estimate  of  the  point-mass 
magnitude.  In  agreement  with  the  previous  classification,  this  adjustment 
entails  four  cases.  The  main  results  are  grouped  in  Table  1,  which  is  self 
explanatory.  This  table  shows  the  error  in  the  tri-diagonal  solution  as  being 
46%  of  the  error  in  the  diagonal  solution.  The  error  in  the  sigma  of  the  tri 
diagonal  solution  amounts  to  about  13%  of  the  error  in  the  sigma  of  the  diagonal 
solution. 

The  error  in  the  five-diagonal  solution  is  much  smaller,  amounting  to  less 
than  23%  of  the  error  in  the  tri  diagonal  solution,  and  thus  to  merely  about  10% 
of  the  error  in  the  diagonal  solution.  With  regard  to  the  error  in  the  sigma  of 
the  five  diagonal  solution,  these  two  numbers  would  be  replaced  by  15%  and  2%, 
respectively.  Although  this  example  is  quite  specialized,  it  nevertheless 
illustrates  that  the  five  diagonal  setup  is  likely  to  match  the  rigorous  setup 


much  more  closely  than  would  be  possible  with  the  tri-diagonal  setup  developed 
in  [Blaha,  1984]. 


Case 

Name 

X 

Error 

a 

X 

Error 

1 

Rigorous 

3.292 

-- 

3.560 

2 

Five  diagonal 

3.049 

-0.243 

3.599 

+  0  039 

3 

Tri -diagonal 

4.367 

+  1.075 

3.821 

+0.261 

4 

Diagonal 

5.622 

+2.330 

1.499 

-2.061 

Table  1 

Least-squares  solutions  and  their  sigmas,  together  with  their  errors, 
in  the  four  weighting  schemes;  the  first  case  is  an  errorless  standard. 


S.  CONCLUSION 


Rigorous  weighting  of  geoldal  observations  along  profiles,  such  as 
encountered,  for  example,  in  satellite  altimetry,  would  be  a  task  exceedingly 
demanding  in  terms  of  computer  run  time.  A  profile  with  m  observation  points 
would  require  an  inversion  of  an  (ix*j  variance  covariance  matrix.  £,  whose 
entries  are  correlated  and  depend  almost  entirely  on  the  geoidal  covariance 
function.  However,  if  the  observations  are  uniformly  distributed,  a  slight 
modification  of  this  function  can  greatly  facilitate  the  formation  of  the  weight 
matrix  of  observations.  P=£  ,  which  subsequently  enters  the  least -squares 
adjustment  process  In  the  suggested  method  this  matrix  Is  five  diagonal .  Its 

six  distinct  entries  are  known  beforehand  for  any  number  of  observations,  which, 
in  addition  to  the  run  time  economies,  also  saves  computer  storage. 

The  pertinent  modification  of  the  covariance  function  is  nearly  negligible, 
especially  if  one  considers  that  this  function  itself  Is  only  approximate,  used 
here  for  a  global  description  of  the  effects  of  the  geoldal  detail  unexpress ible 
by  the  adjustment.  A  five  to  six  times  more  profound  modification  of  the 
covariance  function  (in  terms  of  the  root-mean  square)  leads  to  a  tri-diagonal 
weight  matrix  A  similar  relationship  has  been  noted  in  an  adjustment  of  a 
single  point -mass  magnitude  via  simulated  geoldal  observations,  in  the  sense 
that  the  point  mass  error  associated  with  the  tri-diagonal  weighting  scheme  has 
been  four  to  five  times  larger  than  the  error  associated  with  the  five  diagonal 
weighting  scheme  Since  the  five  diagonal  approach  entails  a  marked  improvement 
in  accuracy  over  its  tri -diagonal  counterpart,  and  since  both  approaches  are 
almost  equally  simple  and  offer  nearly  the  same  computer  economies,  the  former 
Is  preferred  to  the  latter,  and  represents,  in  fact,  the  desired  algorithm  for 
the  adjustment  of  geoidal  profiles  as  sought  in  the  present  study 
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